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1 Introduction 



The study of the bb production cross-section is a powerful test of perturbative quantum chro- 
modynamics (pQCD) calculations. These are available at next-to-leading order (NLO) Q]| and 
with the fixed-order plus next-to-leading logarithms (FONLL) [|2l |3j approximations. In the 
NLO and FONLL calculations, the theoretical predictions have large uncertainties arising from 
the choice of the renormalisation and factorisation scales and the £>-quark mass [4J. Accurate 
measurements provide tests of the validity of the different production models. Recently, the 
LHCb collaboration measured the bb production cross-section in hadron collisions using J/\y 
from b decays O and b — > DfiX decays [6|. The two most recent measurements of the B ± 
production cross-section in hadron collisions have been performed by the CDF collaboration in 
the range pj > 6 GeV/c and |y| < 1 0, where pj is the transverse momentum and y is rapidity, 
and by the CMS collaboration in the range pj > 5 GeV/ c and \y\ < 2.4 [8] . This paper presents 
a measurement of the B ± production cross-section in pp collisions at a centre-of-mass energy 
of yfi = 7 TeV using 34.6 ± 1.2 pb 1 of data collected by the LHCb detector in 2010. The B ± 
mesons are reconstructed exclusively in the B ± — > J/xy^ mode, with J/y — > Both the 

total production cross-section and the differential cross-section, do/dpj, as a function of the 
B ± transverse momentum for < pi < 40 GeV/c and 2.0 < y < 4.5, are measured. 

The LHCb detector [|9l is a single-arm forward spectrometer covering the pseudo-rapidity 
range 2 < rj < 5, designed for the study of particles containing b or c quarks. The detector in- 
cludes a high precision tracking system consisting of a silicon-strip vertex detector surrounding 
the pp interaction region, a large-area silicon-strip detector located upstream of a dipole magnet 
with a bending power of about 4Tm, and three stations of silicon-strip detectors and straw drift- 
tubes placed downstream. The combined tracking system has a momentum resolution Ap/p 
that varies from 0.4% at 5GeV/c to 0.6% at lOOGeV/c, and an impact parameter resolution of 
20 |a.m for tracks with high transverse momentum. Charged hadrons are identified using two 
ring-imaging Cherenkov detectors. Photon, electron and hadron candidates are identified by a 
calorimeter system consisting of scintillating-pad and pre-shower detectors, an electromagnetic 
calorimeter and a hadronic calorimeter. Muons are identified by a muon system composed of 
alternating layers of iron and multiwire proportional chambers. 

The LHCb detector uses a two-level trigger system, the first level (L0) is hardware based, 
and the second level is software based high level trigger (HLT). Here only the triggers used 
in this analysis are described. At the L0 either a single muon candidate with pj larger than 
1.4 GeV/c or a pair of muon candidates, one with pj larger than 0.56 GeV/c and the other with 
pj larger than 0.48 GeV/c, is required. Events passing these requirements are read out and sent 
to an event filter farm for further selection. In the first stage of the HLT, events satisfying one of 
the following three selections are kept: the first one confirms the single-muon candidates from 
L0 and applies a harder pi selection at 1.8 GeV/c; the second one confirms the single-muon 
from L0 and looks for another muon in the event, and the third one confirms the dimuon can- 
didates from L0. Both the second and third selections require the dimuon invariant mass to be 
greater than 2.5 GeV/c 2 . The second stage of the HLT selects events that pass any selections of 
previous stage and contain two muon candidates with an invariant mass within 120 MeV/c 2 of 
the known J/xy mass. To reject high-multiplicity events with a large number of pp interactions, 
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a set of global event cuts (GEC) is applied on the hit multiplicities of sub-detectors. 



2 Event selection 

Candidates for J/\\f — > decay are formed from pairs of particles with opposite charge. 

Both particles are required to have a good track fit quality (j£ 2 /ndf < 4, where ndf represents 
the number of degrees of freedom in the fit), a transverse momentum pj > 0.7 GeV/c and to 
be identified as a muon. In addition, the muon pair is required to originate from a common 
vertex (;£ 2 /ndf < 9). The mass of the reconstructed is required to be in the range 3.04 — 
3.14GeV/c 2 . 

The bachelor kaon candidates used to form B ± — > 7/i//^ ± candidates are required to have 
pj larger than 0.5 GeV/c and to have a good track fit quality (^ 2 /ndf < 4). No particle identi- 
fication is used in the selection of the kaon. A vertex fit is performed that constrains the three 
daughter particles to originate from a common point and the mass of the muon pair to match the 
nominal J/\jr mass. It is required that ^ 2 /ndf < 9 for this fit. To further reduce the combinatorial 
background due to particles produced in the primary pp interaction, only candidates with a de- 
cay time larger than 0.3 ps are accepted. Finally, the fiducial requirement < pj < 40 GeV/c 
and 2.0 < y < 4.5 is applied to the B ± candidates. 



3 Cross-section determination 

The differential production cross-section is measured as 

dg = N B ±(p T ) 

dp T Jzf £totM SS{B ± J/W K± ) ^{J/W -> Ap T ' 

where N B ± (pr) is the number of reconstructed B ± — > J/ylC^ signal events in a given pj bin, Jzf 
is the integrated luminosity, £tot(PT) is me tota l efficiency, including geometrical acceptance, 
reconstruction, selection and trigger effects, 3$(B ± — > J/xj/K^) and 38{Jj\\f — > are the 

branching fractions of the reconstructed decay chain [10], and Apj is the pj bin width. 

Considering that the efficiencies depend on pi and y, we calculate the event yield in bins 
of these variables using an extended unbinned maximum likelihood fit to the invariant mass 
distribution of the reconstructed B ± candidates in the interval 5.15 < M B ± < 5.55 GeV/c 2 . We 
assume that the signal and background shapes only depend on pj. Three components are in- 
cluded in the fit procedure: a Crystal Ball function [|TT| to model the signal, an exponential 
function to model the combinatorial background and a double-Crystal Ball function^] to model 
the Cabibbo suppressed decay B ± — > J/yn . The shape of the latter component is found to fit 
well the distribution of simulated events. The ratio of the number of B ± — > J/xy^ candidates 
to that of the signal is fixed to ^g(B ± 7/v/^ ± ) /^(B ± -> J/yK^ from Ref. [101. The in- 
variant mass distribution of the selected B ± — > 7/i//A' ± candidates and the fit result for one bin 

1 A double-Crystal Ball function has tails on both the low and high mass side of the peak with separate param- 
eters for the two. 
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Figure 1: Invariant mass distribution of the selected B — > J/yrK candidates for one bin 
(5.0 < pj < 5.5 GeV/ c). The result of the fit to the model described in the text is superimposed. 

(5.0 < pr < 5.5 GeV/c) are shown in Fig. [T] The fit returns a mass resolution of 9.14 ± 0.49 
MeV/c 2 , and a mean of 5279.05 ±0.56 MeV/c 2 , where the uncertainties are statistical only. 
Summing over all pj bins, the total number of signal events is about 9100. 

The geometrical acceptance and the reconstruction and selection efficiencies are determined 
using simulated signal events. The simulation is based on Pythia 6.4 generator [T2| with pa- 
rameters configured for LHCb 0"31 . The EvtGen package [fT4l is used to describe the de- 
cays of the B ± and J/y . QED radiative corrections are modelled using PHOTOS 03). The 
Geant4 ffT6l simulation package is used to trace the decay products through the detector. 
Since we select events passing trigger selections that depend on 7/i/a properties only, the trigger 
efficiency is obtained from a trigger-unbiased data sample of J/y events that would still be trig- 
gered if the J/y candidate were removed. The efficiency of GEC is determined from data to be 
(92.6 ±0.3)%, and assumed to be independent of the B ± pi and y. The total trigger efficiency 
is then the product of the J/xy trigger efficiency and the GEC efficiency. The luminosity is mea- 
sured using Van der Meer scans and a beam-gas imaging method 0/7) • The knowledge of the 
absolute luminosity scale is used to calibrate the number of tracks in the vertex detector, which 
is found to be stable throughout the data-taking period and can therefore be used to monitor 
the instantaneous luminosity of the entire data sample. The integrated luminosity of the data 
sample used in this analysis is determined to be 34.6 pb 

The measurement is affected by the systematic uncertainty on the determination of signal 
yields, efficiencies, branching fractions and luminosity. 

The uncertainty on the determination of the signal yields mainly arises from the description 
of final state radiation in the signal fit. The fitted signal yield is corrected by 3.0%, which is 
estimated by comparing the fitted and generated signal yields in the Monte Carlo simulation, 
and an uncertainty of 1.5% is assigned. The uncertainties from the effects of the Cabibbo- 
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suppressed background, multiple candidates and mass fit range are found to be negligible. 

The uncertainties on the efficiencies arise from trigger (0.5 — 6.0% depending on the bin), 
tracking (3.9 — 4.4% depending on the bin), muon identification (2.5%) [5] and the vertex fit 
quality cut (1 .0%). The trigger systematic uncertainty has been evaluated by measuring the trig- 
ger efficiency in the simulation using a trigger-unbiased data sample of simulated 7/i/a events. 
The tracking uncertainty includes two components: the first one is the differences in track recon- 
struction efficiency between data and simulation, estimated with a tag and probe method lfT8ll 
using J/\j/ — > events; the second is due to the 2% uncertainty on the hadronic interaction 

length of the detector used in the simulation. The uncertainties from the effects of GEC, 
mass window cut and inter-bin cross-feed are found to be negligible. The uncertainty due to the 
choice of pi binning is estimated to be smaller than 2.0%. 

The product of 3§{B ± -> J/yfl^) md^(J/y -> jU+jU ) is calculated to be (6.01 ±0.20) x 
10 , by taking their values from Ref. IfTOl with their correlations taken into account. 

The absolute luminosity scale is measured with a 3.5% uncertainty ifTTll . dominated by the 
beam current uncertainty. 

4 Results and conclusion 

The measured B ± differential production cross-section in bins of pj for 2.0 < y < 4.5 is given 
in Table [TJ This result is compared with a FONLL prediction B2l|3l in Fig. |2j A hadronisa- 
tion fraction fi^ B + of (40.1 ± 1.3)% IfTOl is assumed to fix the overall scale of FONLL. The 
uncertainty of the FONLL computation includes the uncertainties on the b-quark mass, renor- 
malisation and factorisation scales, and CTEQ 6.6 f]~9l Parton Density Functions (PDF). Good 
agreement is observed between data and the FONLL prediction. The integrated cross-section is 

oipp^B+X, 0<p T <40 GeV/c, 2.0 < y < 4.5) = 41.4 ± 1.5 (stat.) ±3.1 (syst.) |xb. 

This is the first measurement of B ± production in the forward region at y/s = 7 TeV. 
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Figure 2: Differential production cross-section as a function of the B ± transverse momentum. 
The left plot shows the full pj range, the right plot shows a zoom of the pj range of — 12 
GeV/c. The histogram (left) and the open circles with error bars (right) are the measurements. 
The red dashed lines in both plots are the upper and lower uncertainty limits of the FONLL 
computation. A hadronisation fraction fi^ B + of (40. 1 ± 1 .3) % ifTOl is assumed to fix the overall 
scale. The uncertainty of the FONLL computation includes the uncertainties of the Z?-quark 
mass, renormalisation and factorisation scales, and CTEQ 6.6 PDF. 

Table 1: Differential B production cross-section in bins of pj for 2.0 < y < 4.5. The first and 
second quoted uncertainties are statistical and systematic, respectively. 



PT (GeV/c) 


da/dpT (nb/(GeV/c)) 


Pr (GeV/c) 


do/dp T (nb/(GeV/c)) 


0.0 


-0.5 


1.37 ±0.68 ±0.13 


7.0 


-7.5 


2.42 ±0.20 ±0.18 


0.5 


-1.0 


3.12 ±0.82 ±0.24 


7.5 


-8.0 


2.09 ±0.16 ±0.15 


1.0 


-1.5 


3.90 ± 0.57 ± 0.29 


8.0 


-8.5 


1.44 ±0.11 ±0.11 


1.5 


-2.0 


5.67 ± 1.05 ±0.43 


8.5 


-9.0 


1.33 ±0.11 ±0.10 


2.0 


-2.5 


8.44 ± 1.00 ±0.64 


9.0 


-9.5 


1.22 ±0.10 ±0.09 


2.5 


-3.0 


6.33 ± 0.66 ± 0.48 


9.5- 


-10.0 


0.83 ± 0.08 ± 0.06 


3.0 


-3.5 


5.04 ± 0.45 ± 0.38 


10.0 


-10.5 


0.80 ± 0.08 ± 0.06 


3.5 


-4.0 


6.99 ± 0.68 ± 0.52 


10.5 


-11.0 


0.65 ± 0.07 ± 0.05 


4.0 


-4.5 


5.48 ± 0.47 ± 0.41 


11.0 


-12.0 


0.54 ± 0.04 ± 0.04 


4.5 


-5.0 


6.54 ± 0.79 ± 0.49 


12.0 


-13.0 


0.41 ± 0.04 ± 0.03 


5.0 


-5.5 


4.42 ± 0.44 ± 0.33 


13.0 


-14.5 


0.28 ± 0.02 ± 0.02 


5.5 


-6.0 


4.16 ±0.37 ±0.31 


14.5 


-16.5 


0.17 ±0.02 ±0.01 


6.0 


-6.5 


3.40 ± 0.24 ± 0.25 


16.5 


-21.5 


0.062 ± 0.005 ± 0.005 


6.5 


-7.0 


2.82 ±0.22 ±0.21 


21.5 


-40.0 


0.011 ±0.001 ±0.001 
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